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a View Ten Years after SN 1987A 
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I Abstract. The observation of neutrinos from Supernova 1987A has con- 

firmed the theoretical conjecture that these particles play a crucial role 
during the collapse of the core of a massive star. Only one per cent of 
the energy they carry away from the newly formed neutron star may ac- 
count for all the kinetic and electromagnetic energy responsible for the 
spectacular display of the supernova explosion. However, the neutrinos 
emitted from the collapsed stellar core at the center of the explosion cou- 
^ ■ pie so weakly to the surrounding matter that convective processes behind 

I the supernova shock and/or inside the nascent neutron star might be 

■ required to increase the efficiency of the energy transfer to the stellar 
i mantle and envelope. The conditions for a successful explosion by the 

■ neutrino-heating mechanism and the possible importance of convection 
, in and around the neutron star are shortly reviewed. 
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c/3 . Even after ten exciting years SN 1987A keeps rapidly evolving and develops new, 

' unexpected sides like an aging character. In the first few months the historical 

detection of 24 neutrinos in the underground facilities of the Kamiokande, 1MB, 
and Baksan laboratories caused hectic activity among scientists from very dif- 
ferent fields. In the subsequent years the scene was dominated by the rise and 
I slow decay of light emission in all wavelengths which followed the outbreak of 

the supernova shock and contained a fiood of data about the structure of the 
progenitor star and the dynamics of the explosion. Now that the direct emission 
has settled down to a rather low level, the supernova light which is reflected from 
circumstellar structures provides insight into the progenitor's evolution. Even 
more information about the latter can be expected when the supernova shock 
hits the inner ring in a few years. 

The neutrino detections in connection with SN 1987A were the final proof 
that neutrinos take up the bulk of the energy during stellar core collapse and 
neutron star formation. Lightcurve and spectra of SN 1987A bear clear evidence 
of large-scale mixing in the stellar mantle and envelope and of fast moving Ni 
clumps (see, e.g., J. Spyromilio, D. Wooden, K. Nomoto, this volume). Both 
might indicate that macroscopic anisotropies and inhomogeneities were already 
present near the formation region of Fe group elements during the very early 
stages of the explosion. Spherically symmetric models had been suggesting for 
some time already that regions inside the newly formed neutron star and in the 
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Figure 1. Sketch of the post-collapse stellar core during the neutrino 
heating and shock revival phase. 

neutrino-heated layer around it might be convectively unstable. These theoreti- 
cal results and the observational findings in SN 1987A were motivation to study 
stellar core collapse and supernova explosions with multi-dimensional simula- 
tions. 

In this article convective overturn in the neutrino-heated region around the 
collapsed stellar core is discussed concerning its effects on the neutrino energy 
deposition and its potential importance for the supernova explosion. Convective 
activity inside the proto-neutron star is suggested as a possibly crucial boost of 
the neutrino luminosities on a timescale of a few hundred milliseconds after core 
bounce. The first two-dimensional simulations that follow the evolution of the 
nascent neutron star for more than one second are shortly described. 

2. Neutrino-driven explosions and Convective overturn 

Convective instabilities in the layers adjacent to the nascent neutron star are a 
natural consequence of the negative entropy gradient built up by neutrino heat- 
ing (Bcthe 1990) and are seen in recent two- and three-dimensional simulations 
(Burrows et al. 1995; Herant et al. 1992, 1994; Janka & Miiller 1995, 1996; Mez- 
zacappa et al. 1997; Miller et al. 1993; Shimizu et al. 1994). Although there is 
general agreement about the existence of this unstable region between the radius 
of maximum neutrino heating (which is very close outside the "gain radius" Rg, 
i.e. the radius where neutrino cooling switches into net heating) and the shock 
position Rs, the strength of the convective overturn and its importance for the 
success of the neutrino-heating mechanism in driving the explosion of the star 
is still a matter of vivid debate. 

The effect of convective overturn in the neutrino-heated region on the shock 
is two-fold. On the one hand, heated matter from the region close to the gain 
radius rises outward and at the same time is replaced by cool gas fiowing down 
from the postshock region. Since the production reactions of neutrinos (e^ 
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capture on nucleons and thermal processes) are very temperature sensitive, the 
expansion and cooling of rising plasma reduces the energy loss by reemission 
of neutrinos. Moreover, the net energy deposition by neutrinos is enhanced as 
more cool material is exposed to the large neutrino fluxes just outside the gain 
radius where the neutrino heating rate peaks (the radial dilution of the fluxes 
roughly goes as 1/r^). On the other hand, hot matter floats into the postshock 
region and increases the pressure there. Thus the shock is pushed further out 
which leads to a growth of the gain region and therefore also of the net energy 
transfer from neutrinos to the stellar gas. 

Figure Q displays a sketch of the neutrino cooling and heating regions outside 
the proto-neutron star at the center. The main processes of neutrino energy 
deposition are the charged-current reactions Ue+n — >p+e~ and i^e+P n+e"*". 
The heating rate per nucleon (N) is approximately 
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where Yn and Yp are the number fractions of free neutrons and protons, respec- 
tively, -Li/,52 denotes the luminosity of fe or i/f, in 10^^ erg/s, the radial position 
in 10^ cm, and (£^^5) the average of the squared neutrino energy measured in 
units of 15 MeV. / is the angular dilution factor of the neutrino radiation field 
(the "flux factor" , which is equal to the mean value of the cosine of the angle 
of neutrino propagation relative to the radial direction) which varies between 
about 0.25 at the neutrinosphere and 1 for radially streaming neutrinos far out. 
Using this energy deposition rate, neglecting loss due to re-emission of neutrinos, 
and assuming that the gravitational binding energy of a nucleon in the neutron 
star potential is (roughly) balanced by the sum of internal and nuclear recom- 
bination energies after accretion of the infalling matter through the shock, one 
can estimate the explosion energy to be of the order 
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AM is the heated mass. At the typical heating timescale, Eg^ the (net) total 
gravitational binding energy of the overlying, outward accelerated stellar layers, 
and E'nuc the additional energy from explosive nucleosynthesis which is typically 
a few lO^'^ erg and roughly compensates E'gb for progenitors with main sequence 
masses of less than about 20 M© (see also Burrows, this volume). Since the gain 
radius, shock radius, and At and thus also AM depend on Ly{e^), the sensitivity 
of -Eexp to the neutrino emission parameters is even stronger than suggested by 
Eq. (|. 

In order to get explosions by the delayed neutrino-heating mechanism, cer- 
tain conditions need to be fulfilled. Expansion of the postshock region requires 
sufficiently large pressure gradients near the radius i?cut of the developing mass 
cut. If one neglects self- gravity of the gas in this region and assumes the density 
profile to be a power law, p{r) oc (which is well justified according to numer- 
ical simulations which yield a power law index of n « 3; see also Bethe 1993), 
one gets P{r) oc r^"^^ for the pressure, and outward acceleration is maintained 
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Figure 2. Explosion energies £'>o(i) for ID (dashed) and 2D (solid) 
simulations with different assumed fe and i>e luminosities (labels give 
values in lO^^erg/s) from the proto-neutron star. Below the smallest 
given luminosities the considered 15 Mq star does not explode in ID 
and acquires too low an expansion energy in 2D to unbind the stellar 
mantle and envelope. 



as long as the following condition for the "critical" internal energy density e 
holds: 



GMp/r 
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1 



(n+l)(7-l) 



(3) 



where use was made of the relation P = (7 — 1)5. The numerical value was 
obtained for 7 = 4/3 and n = 3. This condition can be converted into a criterion 
for the entropy per baryon, s. Using the thermodynamical relation for the 
entropy density normalized to the baryon density Ub, s = {e+P)/{ni,T) — J2i Vi^i 
where 77, {i = n, p, e~ , e"^) are the particle chemical potentials divided by the 
temperature, and assuming completely disintegrated nuclei behind the shock 
so that the number fractions of free protons and neutrons are Yp = Yg and 
Y„ = 



Yg, respectively, one gets 
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In this approximate expression a term with a factor Y^ was dropped (its absolute 
value being usually less than 0.5 in the considered region), nucleons are assumed 
to obey Boltzmann statistics, and, normalized to representative values, Mi.i is 
measured in units of 1.1 Mq, pg in 10^ g/cm^, and ry in 10^ cm. Inserting typical 
numbers (T ^ 1.5 MeV, y„ ^ 0.3, i?cut ~ 1-5 • 10^ cm), one obtains s > 15 ks/N, 
which gives an estimate of the entropy in the heating region when the star is 
going to explode. 
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These requirements can be coupled to the neutrino emission of the proto- 
neutron star by the following considerations. A stalled shock is converted into 
a moving one only when the neutrino heating is strong enough to increase the 
pressure behind the shock by a sufficient amount. Considering the Rankine- 
Hugoniot relations at the shock, Bruenn (1993) derived a criterion for the heating 
rate per unit mass, q^, behind the shock that guarantees a positive postshock 
velocity (ni > 0): 

2/3-1 |uo|3 
> /33(/5-l)(7-l) rjR^ ■ 

Here /5 is the ratio of postshock to preshock density, /3 = /Oi/po, 7 the adiabatic 
index of the gas (assumed to be the same in front and behind the shock), and 
T] defines the fraction of the shock radius Rs where net heating by neutrino 
processes occurs: r] = {Rs — Rg)/Rs- uq is the preshock velocity, which is a 
fraction a (analytical and numerical calculations show that typically a ~ l/\/2) 
of the free fall velocity, uq = ay/2GM/r. Assuming a strong shock, one has 
/3= (7 + 1)7(7 — 1) which becomes (3 = 7 for 7 = 4/3. With numbers typical 
of the collapsed core of the 15 star considered by Janka & Miiller (1996), 
i?s = 200km, r] ~ 0.4, and an interior mass M = 1.1 Mq, one finds for the 
threshold luminosities of and i>e'- 

3/2 
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The existence of such a threshold luminosity of the order of 2 • 10^^ erg/s is un- 
derlined by Fig. Q where the explosion energy E^q as function of time is shown 
for numerical calculations of the same post-collapse model but with different 
assumed neutrino luminosities from the proto-neutron star. £'>o is defined to 
include the sum of internal, kinetic, and gravitational energy for all zones where 
this sum is positive (the gravitational binding energies of stellar mantle and 
envelope and additional energy release from nuclear burning are not taken into 
account). For one-dimensional simulations with luminosities below 1.9-10^'^ erg/s 
we could not get explosions when the proto-neutron star was assumed static, and 
the threshold for the z^e and Ue luminosities was 2.2 • lO^'^ erg/s when the neutron 
star was contracting (see Janka &: Miiller 1996). The supporting effects of con- 
vective overturn between the gain radius and the shock described above lead to 
explosions even below the threshold luminosities for the spherically symmetric 
case, to higher values of the explosion energy for the same neutrino luminosi- 
ties, and to a faster development of the explosion. This can clearly be seen by 
comparing the solid (2D) and dashed (ID) lines in Fig. ^. 

The results of Fig. |2| also show that the explosion energy is extremely sensi- 
tive to the neutrino luminosities and mean energies. This holds in ID as well as 
in 2D. Dick McCray in his summary talk of this conference raised the question 
why neutrino-driven explosions should be self-regulated. Which kind of feed- 
back should prevent the explosion from being more energetic than a few times 

gj,g7 Certainly, the neutrino luminosities in current models can hardly power 
an explosion and therefore a way to overpower it is not easy to imagine. Never- 
theless, Fig. ^ and Eq. (^) offer an answer to Dick's question: When the matter 
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Figure 3. Order scheme for the dependence of the post-collapse dy- 
namics on the strength of the neutrino heating as a function of L^{e'^). 



in the neutrino-heated region outside the gain radius has absorbed roughly its 
gravitational binding energy from the neutrino fluxes, it starts to expand out- 
ward (see Eq. (^)) and moves away from the region of strongest heating. Since 
the onset of the explosion shuts off the re-supply of the heating region with cool 
gas, the curves in Fig. ^ approach a saturation level as soon as the expansion 
gains momentum and the density in the heating region decreases. Thus the 
explosion energy depends on the strength of the neutrino heating, which scales 
with the I'e and i^e luminosities and mean energies, and it is limited by the 
amount of matter AM in the heating region and by the duration of the heating 
(see Eq. (^)), both of which decrease when the heating is strong and expansion 
happens fast. 

This also implies that neutrino-driven explosions can be "delayed" (up to a 
few 100 ms after core bounce) but are not "late" (after a few seconds) explosions. 
The density between the gain radius and the shock decreases with time because 
the proto-neutron star contracts and the mass infall onto the collapsed core 
declines steeply with time. Therefore the mass AM in the heated region drops 
rapidly and energetic explosions by the neutrino-heating mechanism become less 
favored at late times. 

Moreover, Fig. |2| tells us that convection is not necessary to get an explosion 
and convective overturn is no guarantee for strong explosions. Therefore one 
must suspect that neutrino-driven type-II explosions should reveal a considerable 
spread in the explosion energies, even for similar progenitor stars. Rotation in 
the stellar core, small differences of the core mass or statistical variations in the 
dynamical events that precede and accompany the explosion may lead to some 
variability. 

The role of convective overturn and its importance for the explosion can be 
further illuminated by considering the three timescales of neutrino heating, Tht, 
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advection of accreted matter through the gain radius into the coohng region and 
onto the neutron star (compare Fig. ||), Tad, and growth of convective overturn, 
Tcv The evolution of the shock — accretion or explosion — is determined by 
the relative sizes of these three timescales. Straightforward considerations show 
that they are of the same order and the destiny of the star is therefore a result 
of a tight competition between the different processes (see Fig. ^ . 

The heating timescale is estimated from the initial entropy Si, the critical 
entropy Sc (Eq. (^)), and the heating rate per nucleon (Eq. (|l|)) as 

^c-^i Sc-gj Rlj{T/2Me\)f 

^ Qt/{kBT) - ^"^^ • 5fcB/iV (L,/2 • 1052erg/s)(e2 ^) ' 



With a postshock velocity ofui = uo/f3^ (7— l)-v/GA//i?s/ (7+I) the advection 
timescale is 

Rs-Rg ( Rg\ -^8,200 

-Tad ~ « 52ms • 1 - -7jr= ' 

ui \ RsJ VMi,i 

where the gain radius can be determined as 

from the requirement that the heating rate, Eq. (|l|), is equal to the cooling rate 
per nucleon, Q' ^ 288(T/2MeV)^ when use is made of the power-law behavior 
of the temperature according to T(r) « 4MeV (Rns/r) with i?ns being the proto- 
neutron star radius (roughly equal to the neutrinosphere radius). The growth 
timescale of convective instabilities in the neutrino-heated region depends on 
the gradients of entropy and lepton number through the growth rate of Ledoux 
convection, ctl {g is the gravitational acceleration): 



dp\ ds f dp\ dY, 



-1/2 

> 50 ms 



.dsjy.^p dr yOYeJs^p dr 

(10) 

The numerical value is representative for those obtained in hydrodynamical sim- 
ulations (e.g., Janka & Miiller 1996). Tcv of Eq. (|lO| ) is sensitive to the detailed 
conditions between neutrinosphere (where Ye has typically a minimum), gain 
radius (where s develops a maximum), and the shock. The neutrino heating 
timescale is shorter for larger values of the neutrino luminosity L^, and mean 
squared neutrino energy (e^), while both r^t and Tad depend strongly on the 
gain radius. Tad also on the shock position. 



3. Convection inside the nascent neutron star 

Convective energy transport inside the newly formed neutron star can increase 
the neutrino luminosities considerably (Burrows 1987). This could be crucial for 
energizing the stalled supernova shock (Mayle & Wilson 1988; Wilson & Mayle 
1988, 1993; see also Sect. |]). 
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Figure 4. Left: Luminosities Li,{f?) and mean energies (ejy)(i) of i/g 
and for a 1.1 proto-neutron star without ("ID"; dotted) and 
with convection ("2D"; solid). Time is measured from core bounce. 
Right: Angular variations of the neutrino flux at different times for the 
2D simulation. 



Recent two-dimensional simulations by Keil et al. (1996, 1997) and Keil 
(1997) have followed the evolution of the proto-neutron star formed in the core 
collapse of a 15 star for a period of more than 1.2 seconds. The simulations 
were performed with the hydrodynamics code Prometheus. A general relativistic 
ID gravitational potential with Newtonian corrections for asphericities was used, 
$ = <&fj3^ + (*&^D — ^I^id), and a flux-limited (equilibrium) neutrino diffusion 
scheme was applied for each angular bin separately ("l^D"). 

The simulations show that convectively unstable surface-near regions (i.e., 
around the neutrinosphere and below an initial density of about 10'''^g/cm^) 
exist only for a short period of a few ten milliseconds after bounce, in agreement 
with the findings of Bruenn & Mezzacappa (1994), Bruenn et al. (1995), and 
Mezzacappa et al. (1997). Due to a flat entropy profile and a negative lepton 
number gradient, convection, however, also starts in a layer deeper inside the 
star, between an enclosed mass of 0.7 and 0.9 Mq, at densities above several 
lO^^g/cm'^. From there the convective region digs into the star and reaches 
the center after about one second. Convective velocities as high as 5 • lO^cm/s 
are reached (about 10-20% of the local sound speed), corresponding to kinetic 
energies of up to 1—2 • 10^*^ erg. Because of these high velocities and rather flat 
entropy and composition profiles in the star, the overshoot (and undershoot) 
regions are large. 

The coherence lengths of convective structures are of the order of 20-40 
degrees (in 2D!) and coherence times are of the order of 10 ms which corre- 
sponds to only one or two overturns. The convective pattern is therefore very 
time-dependent and nonstationary. Convective motions lead to considerable 
variations of the composition. The lepton fraction (and thus the abundance of 
protons) shows relative fluctuations of several 10%. The entropy differences in 



8 



rising and sinking convective bubbles are much smaller, only a few per cent, 
while temperature and density fluctuations are typically less than one per cent. 

The energy transport in the neutron star is dominated by neutrino diffu- 
sion near the center, whereas convective transport plays the major role in a 
thick intermediate layer where the convective activity is strongest, and radiative 
transport takes over again when the neutrino mean free path becomes large near 
the surface of the star. But even in the convective layer the convective energy 
flux is only a few times larger than the diffusive flux. This means that neutrino 
diffusion can never be neglected. 

There is an important consequence of this latter statement. The convective 
activity in the neutron star cannot be described and explained as Ledoux convec- 
tion. Applying the Ledoux criterion for local instability, Ci^{r) = {p/g)ai^ > 
with (Tl from Eq. (10) and YJ= replaced by the total lepton fraction Yiep in the 
neutrino-opaque interior of the neutron star, one finds that the convecting re- 
gion should actually be stable, despite of slightly negative entropy and lepton 
number gradients. In fact, below a critical value of the lepton fraction (e.g., 
^iep,c = 0.148 for p = 10^^ g/cm^ and T = 10.7 MeV) the thermodynamical 
derivative {dp/dYiep)s^p changes sign and becomes positive because of nuclear 
and Coulomb forces in the high-density equation of state. Therefore negative 
lepton number gradients should stabilize against convection in this regime. How- 
ever, an idealized assumption of Ledoux convection is not fulfilled in the situ- 
ations considered here: Because of neutrino diffusion energy exchange and, in 
particular, lepton number exchange between convective elements and their sur- 
roundings are not negligible. Taking the neutrino transport effects on liep into 
account in a modified Quasi-Ledoux criterion (Keil 1997 and Keil et al. 1997) 
one predicts instability exactly where convective action happens in the two- 
dimensional simulation. 



4. Conclusions 

Convection inside the proto-neutron star can raise the neutrino luminosities 
within a few hundred ms after core bounce (Fig. ^). In the considered collapsed 
core of a 15 Mq star L^^ and Lp^ increase by up to 50% and the mean neutrino 
energies by about 15% at times later than 200-300 ms post bounce. This favors 
neutrino-driven explosions on timescales of a few hundred milliseconds after 
shock formation. Also, the deleptonization of the nascent neutron star is strongly 
accelerated, raising the Ve luminosities relative to the I'e luminosities during this 
time. This helps to increase the electron fraction Yg in the neutrino-heated 
ejecta and might solve the overproduction problem of = 50 nuclei during 
the early epochs of the explosion (Keil et al. 1996). Anisotropic mass motions 
due to convection in the neutron star lead to gravitational wave emission and 
anisotropic radiation of neutrinos. The angular variations of the neutrino flux 
determined by the 2D simulations are of the order of 5-10% (Fig. ^). With the 
typical size and short coherence times of the convective structures, however, the 
global anisotropy of the neutrino emission from the cooling proto-neutron star 
is certainly less than 1% (more likely only 0.1%, since in 3D the structures tend 
to be smaller) and kick velocities in excess of 300 km/s can definitely not be 
explained. 
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In more recent simulations, mass accretion and rotation of the forming 
neutron star were included. Rotation has very interesting consequences, e.g., 
leads to a suppression of convective motions near the rotation axis because of 
a stabilizing stratification of the specific angular momentum, an effect which 
can be understood by applying the Solberg-H0iland criterion for instabilities 
in rotating, self-gravitating objects. Future simulations will have to clarify the 
influence of the nuclear equation of state on the presence of convection in nascent 
neutron stars. Also, a more accurate treatment of the neutrino transport in 
combination with a statc-of-thc-art description of the neutrino opacities of the 
nuclear medium is needed to confirm the existence of a convective episode during 
neutron star formation and to study its importance for the explosion mechanism 
of type-II supernovae. 
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